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PREFACE 


This  study  was  performed  by  Mr.  R.  G.  Willey  with  the  technical  assistance 
of  Mr.  Keith  Knight.  Mr.  Don  Smith  of  Resource  Management  Associates  provided 
advice  during  critical  parts  of  the  study.  The  study  was  managed  under  the 
direction  of  Dr.  Richard  Punnett  of  the  Huntington  District  who  was  also 
responsible  for  providing  all  of  the  data  necessary  for  calibrating  and 
verifying  the  UEC-  o  computer  model. 

The  HEC-5Q  computer  model  was  originally  developed  under  funding  from  the 
Corps'  Environmental  and  Water  Quality  Operational  Studies  (EWQOS)  program. 

This  study  serves  as  a  field  verification  demonstration  of  HKC-5Q. 

Mr.  Vern  Bonner  is  Chief  of  the  Training  Division  and  Mr.  Willey's 
supervisor.  Mr.  Bill  Eichert  is  Director  of  the  Hydrologic  Engineering  Center. 
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1.  INTRODUCTION 


1.1  PURPOSE 

~  >>The  primary  purpose  of  this  study  is  to  apply  newly  developed  computer 
technology  to  water  quality  analysis  of  the  reservoir  system  network  in  the 
Kanawha  River  Basin,  West  Virginia.  In  Hay  1985,  the  Huntington  District 
initiated  a  technical  assistance  project  with  the  Hydrologic  Engineering 
Center  (HEC)  to  apply  HEC-5Q  to  the  Kanawha  River  Basin;  this  report 
summarizes  the  results  of  the  project  and  is  intended  as  a  reference  for  the 
Huntington  District. 

A  secondary  purpose  of  this  report  is  to  provide  an  example  that  will 
serve  as  a  guide  for  those  who  are  applying  the  HEC-5Q  computer  program  to 
other  river  basins.  Although  the  details  of  the  applications  will  change,  the 
process  described  is  expected  to  remain  the  same. 


1 . 2  SCOPE 


'■1  In  1979,  the  HEC  began  to  develop  the  HEC-5Q  computer  program  to  provide 
better  system  water  quality  analysis  capabilities  in  support  of  the  Corps* 
water  control  management  program.  The  focus  of  this  program  is  to  evaluate 
water  quality  conditions  associated  with  proposed  plans  of  reservoir  system 
regulation  or  to  determine  improved  reservoir  system  operations  for  multiple 
project  purposes  including  water  quality.  T  vs.' 1 1 

The  computer  program  HEC-5Q  is  described  in  the  Users  Manual  (HEC,  1984a] 
in  addition  to  several  technical  papers  {Duke,  1984  and  Willey:  1983,  1985a 
and  1986]  and  reports  [Willey,  1985b].  ^ 

This  report  identifies  the  type  of  data  necessary  to  execute  HEC-5Q,  where 
it  is  typically  found  and  how  it  is  manipulated  prior  to  input  to  the 
program.  It  also  describes  the  calibration  and  verification  process,  and 
finally  shows  the  results  of  the  executions. 

1 . 3  AUTHORITY 

The  study  was  authorized  by  the  Huntington  District,  with  DA  Form  2544 
No.  E86-85-HW-33  dated  29  May  1985.  An  increase  in  funds  was  authorized  with 
a  DA  Form  2544,  Change  Order  Number  1,  dated  11  October  1985.  Hatching  funds 
were  provided  by  the  Office  of  the  Chief  of  Engineers  as  a  demonstration 
project  for  Corps-wide  technology  transfer. 

1 . 4  PROCEDURE 

The  study  was  initiated  with  a  conference  between  staff  from  Huntington 
District,  Ohio  River  Division,  Office  of  the  Chief  of  Engineers  and  the 
Hydrologic  Engineering  Center  (HEC).  A  study  proposal  prepared  by  the  HEC  was 
discussed  and  a  briefing  on  the  HEC-5Q  model  was  presented.  The  District 
decided  to  accept  the  proposal  and  have  HEC  begin  work  in  June  1985. 


The  first  step  of  the  procedure  involved  data  assembly.  The  District 
personnel  were  the  most  knowledgeable  about  data  availability  and  were 
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therefore  responsible  for  obtaining  the  required  data.  The  HI£C,  being  more 
familiar  with  the  model  input,  assumed  responsibility  for  requesting  the  data 
necessary  to  execute  the  model. 


The  second  step  involved  the  manipulation  of  some  data  by  utility  programs 
prior  to  use  as  input  to  HEC-SQ.  Usually  meteorological  data  from  the 
National  Climatic  Center  and  channel  cross-section  data  are  processed  through 
independent  utility  programs  (i.e.,  WEATHER,  HEATX,  and  GEOA)  which  structure 
the  data  in  the  proper  units  and  format  required  by  HEC-5Q.  This  task  was 
performed  by  HEC. 

The  third  task  involved  model  calibration.  Model  calibration  is  the  most 
difficult  step  in  the  procedure.  It  involves  modification  of  model  inputs  for 
each  trial  execution.  The  model  inputs  which  usually  change  are  parameters 
which  are  not  easily  measurable  nor  available.  Examples  of  these  include  the 
reservoir  diffusion  coefficients,  the  distribution  of  absorbed  radiation,  and 
the  estimated  variation  both  spatially  and  temporally  of  the  inflowing  water 
quality  concentrations  at  all  boundaries.  This  task  was  performed  by  HEC  with 
expert  technical  assistance  by  Mr.  Donald  J.  Smith,  a  consultant  engineer  with 
Resource  Management  Associates  of  Lafayette,  California. 

Step  four  was  model  verification.  Model  verification  involved  applying 
the  model  to  an  independent  period  of  data  without  modification  of  the 
previously  mentioned  calibration  parameters. 

All  four  steps  in  the  procedure  are  discussed  in  detail  in  Section  4  of 
this  report. 
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2 .  BASIS  DESCRIPTION* 


2.1  PHYSIOGRAPHY 

The  Kanawha  River  basin  drains  12,300  square  miles  in  North  Carolina, 
Virginia,  and  West  Virginia.  It  is  a  tributary  to  the  Ohio  River  at  Pt. 
Pleasant,  West  Virginia,  as  shown  in  Figure  2.1.  The  headwaters  of  the  basin 
are  in  the  northwest  comer  of  North  Carolina.  From  there,  the  New  River 
flows  about  340  miles  through  the  southwest  comer  of  Virginia  and  into 
southcentral  West  Virginia.  Drainage  from  the  Greenbrier  River,  which  flows 
along  the  eastern  border  of  West  Virginia,  joins  the  New  River  near  Hinton. 

The  Gauley  River,  which  drains  central  West  Virginia,  joins  the  New  River  near 
Kanawha  Falls.  The  confluence  of  the  New  and  Gauley  Rivers  forms  the  Kanawha 
River  which  then  flows  97  miles  to  the  Ohio  River.  The  Kanawha  River  is 
joined  by  the  Elk  River  at  Charleston,  which  also  drains  central  West 
Virginia.  The  drainage  areas  of  these  tributaries  are  listed  in  Table  2.1. 

The  area  of  concern  in  this  report  includes  the  Elk  and  Gauley  River  drainage 
and  the  New  River  from  Hinton  downstream  to  Winfield  Lock  and  Dam  (23  miles 
downstream  of  Charleston)  on  the  Kanawha  River. 

The  channels  of  the  New  River  and  its  tributaries  are  confined  by 
mountainous  terrain  to  narrow  valleys.  As  a  result,  the  flood  plain  is  often 
only  slightly  wider  than  the  low-water  channel  which  varies  in  width  from  200 
to  1,000  feet.  In  many  locations,  banks  are  sheer  bluffs  rising  from  the  edge 
of  the  river.  The  channels  of  the  Gauley  and  Elk  Rivers  are  also 
characterized  by  narrow  valleys  and  gorges.  The  main  channel  from  Kanawha 
Falls  to  Winfield  consists  of  lock-and-dam  navigation  pools.  In  this  reach, 
banks  of  the  river  are  relatively  high,  ranging  from  25  to  55  feet  above  low 
water.  The  width  of  the  river  at  normal  stage  varies  from  500  to  800  feet. 

Topography  of  the  Kanawha  basin  is  highly  dissected  with  local  relief 
ranging  from  about  350  to  1500  feet.  Surface  slopes  in  excess  of  25%  are 
common.  The  elevation  of  the  ridge  bordering  the  watershed  ranges  from  about 
4000  feet  in  the  south,  to  2000  feet  in  the  east,  to  800  feet  in  the  northwest 
near  the  Ohio  River.  Most  of  the  mountains  are  forested. 

Table  2.1 


Kanawha  River  Basin 
Principal  Rivers 


River 

Drainage 

Area 

(souare  miles) 

Elevation 
of  Source 
(feet) 

Elevation 
of  Mouth 
(feet) 

Average 

Slope 

(feet/mile) 

New 

6920 

3800 

652 

9.2 

Kanawha 

12300 

652 

538 

1.2 

Greenbrier 

1653 

26  75 

1360 

87.4 

Gauley 

440 

4000 

652 

31.0 

Elk 

1536 

3600 

563 

16.8 

*Partially  from  Real-Time  Flood  Forecasting  and  Reservoir  Control  for  the 
Kanawha  River  Basin  [HEC,  1984b] 
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POINT  PLEASANT 


FIGURE  2.1 

KANAWHA  RIVER  BASIN 


As  8ho%m  in  Table  2.1,  streambed  slopes  of  the  principal  rivers  vary  from 

1.2  to  87.4  feet  per  mile.  Slopes  of  the  tributaries  are  much  higher;  most 
are  in  excess  of  100  feet  per  mile  and  many  exceed  400  feet  per  mile. 

2.2  CLIMATOLOGY 

The  Kanawha  River  Basin  lies  beyond  the  immediate  climatic  effect  of  the 
Atlantic  Ocean,  and  is  therefore  characterized  by  a  marked  temperature 
contrast  between  summer  and  winter.  The  basin  is  affected  by  prevailing 
westerly  winds  which  are  frequently  interrupted  by  northward  and  southward 
surges  of  relatively  warm  and  cold  air,  respectively.  These  atmospheric 
movements  are  accompanied  by  the  passage  of  high-  and  low-pressure  areas  which 
move  in  a  generally  easterly  direction  across  the  United  States  in  the  colder 
half  of  the  year.  From  May  through  October,  the  basin  is  affected  by  showers 
and  thunderstorms  that  occur  in  the  broad  current  of  air  that  tends  to  sweep 
northeastward  from  the  Gulf  of  Mexico.  Winter  climate  throughout  the  basin  is 
moderately  severe  with  frequent  alternations  of  fair  and  stormy  weather. 

Summer  is  marked  by  hot  and  showery  weather,  with  cooler  temperatures  and 
frequent  thunderstorms  in  the  mountains.  Thunderstorms  occur  on  an  average  of 
40  to  50  days  per  year,  with  the  greatest  frequency  in  June  and  July. 


2.3  RESERVOIR  SYSTEM 

Three  of  the  four  major  storage  impoundments  in  the  Kanawha  River  Basin 
are  operated  by  the  Corps  of  Engineers.  The  Corps  projects  are  Sutton  Lake  on 
the  Elk  River,  Summersville  Lake  on  the  Gauley  River,  and  Bluestone  Lake  on 
the  Mew  River.  The  fourth  project  is  Clay tor  Dam,  a  hydropower  project  owned 
and  operated  by  Appalachian  Power  Company.  Clay tor  Dam  is  above  Bluestone 
Lake  Dam  and  therefore  of  limited  interest  in  this  study.  Basic 
characteristics  of  the  reservoirs  are  presented  in  Table  2.2. 

Table  2.2 

Kanawha  River  Basin 
Major  Reservoirs 


Flood  Control  Storage 


Reservoir 

Drainage  Area 

Winter* 

Summer* 

Name 

(square  miles) 

(ac-ft) 

(ac-ft) 

Claytor  Reservoir 

2,382 

** 

** 

Bluestone  Lake 

4,603 

600,100 

592,600 

Summersville  Lake 

803 

355,524 

221,884 

Sutton  Lake 

543 

236,130 

205,570 

*  Different  winter  and  summer  storages  reflect  the  seasonal  potential  for 
runoff.  Winter  reservation  is  for  the  period  1  December  through  15  April. 

**  An  active  power  storage  of  95,000  ac-ft  (.75  inches)  provides  incidental 
flood  control  during  minor  floods. 


8 


The  Kanawha  River  Basin  reservoirs  are  operated  to  control  flows  at 
designated  control  points  on  the  Kanawha  and  its  tributaries.  Control  points 
were  selected  based  on  flood-damage  surveys,  flood  operation  procedures,  and 
availability  of  reliable  streamflow  rating  stations  to  permit  satisfactory 
data  collection.  The  control  points  and  their  pertinent  characteristics  are 
shown  in  Table  2.3.  A  schematic  of  the  Kanawha  River  Basin  reservoirs  and 
control  points  are  shown  in  Figure  2.2. 


Table  2.3 


Kanaifha 

River  Basin 

Control  Points 

Control  Point 

Channel 

Capacity  (cfs) 

Location 

Stream 

Winter* 

Summer* 

Sutton  Dam 

Elk  River 

12,000 

8,000 

Frametown 

Elk  River 

17,600 

14,800 

Clay 

Elk  River 

35,200 

35,200 

Queen  Shoals 

Elk  River 

28,600 

28,600 

Summersville  Dam 

Gauley  River 

20,500 

20,500 

Belva 

Gauley  River 

42,500 

42,500 

Hinton 

Mew  River 

90.200 

90,200 

Kanawha  Falls 

Kanawha  River 

146,000 

146,000 

London  L/D 

Kanawha  River 

150,000 

150,000 

Marmet  L/D 

Kanawha  River 

150,000 

150,000 

Charleston** 

Kanawha  River 

150,000 

150,000 

Winfield  L/D 

Kanawha  River 

45*** 

45*< 

*  Different  winter  and  summer  capacities  reflect  the  greater  potential  for 
damages  to  agricultural  areas  during  the  growing  season.  Winter  is  from  1 
December  through  15  April. 

**  Old  Lock  So.  6  gage  with  a  slope  rating;  control  discharge  estimated. 

***  Stage  in  feet. 
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Figure  2.2 

KANAWHA  RIVER  BASIN  SCHEMATIC  * 


*From  HEC-5Q:  System  Water  Qaulity  Ilodeling  [Willey,  1986] 
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HEC-5Q  RESERVOIR  REGULATION  MODEL 


The  "HEC-5Q,  Simulation  of  Flood  Control  and  Conservation  Systems 
(Including  Water  Quality  Analysis)*'  computer  model  [HEC,  1984a]  has  been 
developed  specifically  for  evaluating  the  type  of  problem  shown  in  Figure 
2.2.  The  model  is  capable  of  evaluating  a  reservoir  system  of  up  to  ten 
reservoirs  and  thirty  control  points.  The  model  will  analyze  a  best  system 
operation  for  water  quantity  and  quality  by  evaluating  operational  concerns 
like  flood  control,  hydropower,  water  supply,  irrigation  diversions  and  other 
instream  flow  needs. 

3.1  FLOW  SIMULATION  MODULE 

In  addition  to  assisting  in  the  sizing  of  the  flood  control  and 
conservation  storage  requirements  for  each  recommended  project,  the  flow 
simulation  module  was  developed  to  assist  in  planning  studies  for  evaluating 
proposed  reservoirs  in  a  system.  The  program  can  also  be  useful  for  selecting 
proper  reservoir  operational  releases  in  order  to  maximize  water  quality 
conditions  or  other  operational  downstream  targets. 

3.2  WATER  QUALITY  SIMULATION  MODULE 

The  water  quality  simulation  module  was  developed  to  simulate  temperature, 
as  well  as  three  user-selected  conservative  and  three  user-selected 
non-conservative  constituents.  The  model  also  allows  dissolved  oxygen  to  be 
simulated  if  the  user  selects  either  carbonaceous  or  nitrogenous  oxygen 
demanding  constituents.  An  option  for  phytoplankton  evaluation  is  also 
available. 

The  water  quality  simulation  module  accepts  system  flows  generated  by  the 
flow  simulation  module  and  computes  the  distribution  of  all  the  water  quality 
constituents  in  the  reservoirs  and  their  associated  downstream  reaches.  The 
program  also  selects  the  gate  openings  for  reservoir  withdrawal  structures  to 
meet  user-specified  water  quality  objectives  at  downstream  control  points. 

With  these  capabilities,  the  planner  may  evaluate  the  effects  of  proposed 
reservoir-stream  modifications  on  water  quality  and  determine  how  a  reservoir 
intake  structure  should  be  operated  to  achieve  desired  water  quality 
objectives  within  the  system. 

The  reservoirs  are  represented  conceptually  by  one-dimensional  horizontal 
elements.  Each  horizontal  element  is  characterized  by  an  area,  thickness  and 
volume.  In  the  aggregate,  the  assemblage  of  layered  volume  elements  is  a 
geometric  representation  of  the  prototype  reservoir.  This  one-dimensional 
representation  has  been  sho%m  to  adequately  represent  water  quality  conditions 
in  many  deep,  well-stratified  reservoirs  by  Baca  [1977]  and  Water  Resources 
Engineers  [1968,  1969a,  1969b]. 

Each  horizontal  layer  is  assumed  to  be  completely  mixed  with  all  isopleths 
parallel  to  the  water  surface  both  laterally  and  longitudinally.  External 
inflows  and  withdrawals  occur  as  sources  or  sinks  within  each  layer  and  are 
instantaneously  dispersed  and  homogeneously  mixed  throughout  each  element  from 
the  headwaters  of  the  impoundment  to  the  dam.  It  is  not  possible,  therefore, 
to  look  at  longitudinal  variations  in  water  quality  constituents.  Simulation 
results  are  most  representative  of  conditions  in  the  main  reservoir  body. 
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Vertical  advection  is  governed  by  the  location  of  inflow  to,  and  outflow 
from,  the  reservoir.  Thus,  the  computation  of  the  zones  of  distribution  and 
withdrawal  for  inflows  and  outflows  are  of  considerable  significance  in 
operation  of  the  model.  The  WES  withdrawal  method  [Bohan,  1973]  is  used  for 
determining  the  allocation  of  outflow.  The  Debler  inflow  allocation  method 
[Debler,  1959]  is  used  for  the  placement  of  inflows. 

Vertical  advection  is  the  net  interelement  flow  and  is  one  of  two 
transport  mechanisms  used  in  the  module  to  transport  water  quality 
constituents  between  elements.  Effective  diffusion  is  the  other  transport 
mechanism.  The  effective  diffusion  is  composed  of  molecular  and  turbulent 
diffusion  and  convective  mixing. 

The  stream  system  is  represented  conceptually  as  a  linear  network  of 
segments  or  volume  elements.  Each  element  is  characterized  by  length,  width, 
cross  sectional  area,  hydraulic  radius.  Manning's  n  and  a  flow  and  depth 
relationship.  Flow  rates  at  stream  control  points  are  calculated  within  the 
flow  simulation  module  using  any  one  of  several  programmed  hydrologic  routing 
methods.  Within  the  flow  simulation  module,  incremental  local  flows  (i.e., 
inflow  between  adjacent  control  points)  are  assumed  to  be  located  at  the 
nearest  control  point. 

Within  the  water  quality  simulation  module,  the  incremental  local  flow  may 
be  divided  into  components  and  placed  at  different  locations  within  the  stream 
reach  (i.e.,  that  portion  of  the  stream  bounded  by  the  two  control  points).  A 
flow  balance  is  used  to  determine  the  flow  rate  at  element  boundaries.  Any 
flow  imbalance  (i.e.,  the  difference  between  the  flow  at  the  upstream  control 
point,  including  all  tributary  inflows,  and  the  flow  at  the  doiimstream  control 
point)  is  distributed  uniformly  to  the  flows  at  each  element  boundary.  Once 
interelement  flows  are  established,  the  depth,  surface  width,  and  cross 
sectional  area  are  computed  at  each  element  boundary  assuming  normal  depth. 

3 . 3  GATE  SELECTION 

Once  the  desired  reservoir  release  and  the  target  water  quality  needed  to 
meet  downstream  requirements  have  been  computed,  the  gate  selection  algorithm 
determines  which  ports  should  be  open  and  what  flow  rate  should  pass  through 
each  open  port  in  order  to  maximize  a  particular  function  of  the  downstream 
water  quality  target  concentrations.  Solution  of  this  problem  is  accomplished 
by  using  mathematical  optimization  techniques.  The  objective  function  is 
related  to  meeting  downstream  target  qualities  subject  to  various  hydraulic 
constraints  on  the  individual  ports. 

The  algorithm  considers  a  sequence  of  problems,  each  representing  a 
different  combination  of  open  ports.  For  each  combination,  the  optimal 
allocation  of  total  flow  to  ports  is  first  determined;  and  then  a  water 
quality  index  is  determined  for  the  optimal  allocation  of  flows.  The 
combination  of  open  ports  with  the  highest  water  quality  index  and  its 
associated  allocation  of  flows  defines  the  optimal  operation  strategy  for  the 
time  period  under  consideration. 
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3.4  FLOW  ALTBSATIONS 


The  flow  alteration  routine  is  designed  to  change  the  reservoir  releases, 
computed  by  the  flow  simulation  module,  to  better  satisfy  the  stream  control 
point  water  quality  objectives.  The  routine  is  designed  about  a  mass  balance 
for  all  reservoir  releases  and  all  control  points  affected  by  those  releases. 
Tributary  inflows  and  other  flow  changes  are  included.  Second  order  effects, 
such  as  changes  in  reaeration  and  external  heating  due  to  increased  or 
decreased  stream  surface  area,  are  not  included. 

Thus  the  flow  augmentation  requirement  can  be  computed  for  each  control 
point  and  for  each  constituent.  The  various  computed  flow  rates  are  then 
combined  by  using  the  coefficients  of  the  linear  programming  objective 
function  and  the  deviation  of  the  respective  constituent  concentrations  from 
the  target  concentrations  at  each  respective  control  point. 

Once  the  flow  augmentation  requirement  is  determined,  the  flow  simulation 
module  is  recalled;  and  the  daily  computations  for  flow  and  water  quality  are 
solved  again  for  the  final  results.  On  this  simulation,  second  order  effects 
are  included  in  the  analysis. 

3.5  SUMMARY 

The  HEC-5Q  model  is  capable  of  simulating  the  operational  effects  of  as 
many  as  ten  reservoirs  on  the  stream  network  of  the  basin.  Each  reservoir  may 
be  operated  to  satisfy  a  number  of  objectives  including  flood  control,  low 
flow,  hydropower  production,  water  supply  and  water  quality  control.  The 
water  quality  portion  of  the  model  will  simulate  temperature  and  eight  water 
quality  constituents  including  dissolved  oxygen  and  phytoplankton.  The  model 
will  internally  determine  the  water  quality  needed  from  all  reservoir  releases 
to  meet  specified  downstream  water  quality  objectives  and  will  determine  the 
gate  openings  in  each  reservoir  that  will  yield  the  appropriate  reservoir 
release  water  quality.  Should  it  be  necessary,  flows  will  be  altered  to 
ensure  that  downstream  water  quality  objectives  are  met.  The  model  selects 
the  best  solution  for  system-wide  reservoir  operation  on  a  daily  basis. 
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4.  METHOD  OF  STUDY 


Input  data  to  simulate  reservoir  regulation  for  evaluation  of  water 
quality  impacts  was  developed  based  on  the  generalized  reservoir  system 
simulation  program  HEC-5Q  [HEC,  1984al 

Development  of  the  HEC-5Q  input  data  comprised  the  following  tasks: 

1)  Determining  the  basin  network  configuration  needed  to  adequately 
describe  the  reservoir  system,  as  well  as  an  appropriate  simulation 
period  and  time  interval; 

2)  Cotopiling  project  data  describing  the  reservoir  system; 

3}  Selecting  low-flow  events  for  simulation;  compiling  time  series  data 
on  the  selected  low-flow  events;  and  calibrating  the  data  models  by 
successive  trial  and  error  adjustments  of  appropriate  variables. 

4.1  BASIN  NETWORK  AND  TIME  CONSIDERATIONS 

The  basin  network  was  identical  in  design  to  the  network  previously  used 
in  the  Real-Time  Flood  Forecasting  Study  report  [HEC,  1984b]  except  that  the 
upstream  boundary  on  the  New  River  was  Hinton  and  the  downstream  boundary  on 
the  Kanawha  River  was  Winfield  Lock  and  Dam.  The  Elk  and  Gauley  River 
analysis  began  with  knovm  inflow  to  Sutton  and  Summersville  Lakes 
respectively.  The  schematic  of  the  network  is  shown  in  Figure  2.2.  The 
reason  for  different  boundary  limits  on  the  New/Kanawha  Rivers  can  be 
attributed  to  a  District  decision  based  on  interests  in  water  quality 
conditions . 

The  HEC-5Q  model  presently  simulates  either  daily  or  monthly  data.  This 
study  used  the  daily  time  interval.  The  District  selected  four  years  of 
specific  water  quality  interest:  1976,  1979,  1980  and  1983.  The  season  of 
interest  was  June  through  September. 

The  data  for  1976  was  unavailable  prior  to  October,  so  a  trial  analysis 
was  performed  for  October  1976.  This  case  had  a  minimum  of  calibration  data 
available. 

The  data  for  1979  and  1980  were  sufficient  for  viable  calibration 
atten^ts.  The  data  for  1983  was  not  readily  available  for  calibration  but  was 
made  available  at  a  later  date  for  verification  analysis. 

4.2  PROJECT  DATA  FOR  RESERVOIR  SYSTEM 

Project  data  defines  the  pertinent  physical  and  operational 
characteristics  of  the  prototype  reservoir  system  which  are  relati-vely  fixed 
in  time.  For  example,  project  data  include  reservoir  storage,  surface  area, 
and  pool  elevation  relationships;  outflow  and  elevation  relationships;  and 
routing  par2uneter8.  Certain  operation  parameters  are  also  relatively  static, 
such  as  storage  allocations  and  channel  capacities. 
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The  primary  sources  of  project  data  were  the  reservoir  regulation  manuals 
for  the  individual  reservoirs  [U.S.  Army,  1953,  1961,  1966a]  and  the  master 
manual  for  reservoir  regulation  [U.S.  Army,  1966b]. 

4.2.1  Storage,  Elevation  and  Outflow 

Relationships  between  storage,  area,  outflow  and  elevation  are  represented 
in  HEC-5Q  by  tables  of  corresponding  linearly  interpolated  points.  Points 
were  selected  to  obtain  acceptable  piecewise  linear  functions  to  represent 
curvilinear  relationships.  Outflow  capabilities  were  specified  as  the  sum  of 
spillway  and  main  outlet  capacities  where  appropriate.  Sutton,  which  has  a 
gated  spillway,  has  full  combined  sluice  gate  and  spillway  capacity  for 
elevations  above  the  spillway  crest.  Summersville,  which  has  an  overflow 
spillway,  has  total  main  outlet  capacity  below  the  spillway  crest;  for 
elevations  above  the  spillway  crest,  the  outflow  capability  is  specified  as 
the  spillway  plus  any  available  main  outlet  capability. 

4.2.2  Storage  Allocations 

Reservoir  storage  allocations  were  assigned  among  five  zones:  inactive, 
buffer,  conservation,  flood  control,  and  emergency  flood  surcharge  storage. 
Zones  were  defined  by  storage  values  at  the  tops  of  the  zones.  Seasonal  pool 
variations  of  conservation  pools  were  included  where  appropriate.  Some  of  the 
buffer  or  inactive  storage  levels  were  not  specifically  defined  in  project 
manuals;  consequently,  values  were  estimated. 

4.2.3  Control  Points  and  Channel  Capacities 

The  control  points  and  their  channel  capacities  are  shown  in  Table  2.3. 

The  channel  capacities  were  obtained  from  the  reservoir  regulation  manuals  for 
the  individual  reservoirs,  when  available.  At  Summersville,  the  channel 
capacity  was  made  equal  to  the  maximum  capacity  of  the  three  main  outlet 
valves,  about  20,500  cfs,  as  specified  in  the  manual.  A  seasonal  channel 
capacity  was  specified  for  Sutton  in  accordance  with  its  manual. 

4.2.4  Minimum  Flows 

Minimum  releases  specified  in  project  manuals  for  the  Corps  reservoirs 
were  included  in  the  HEC-5Q  model.  The  minimum  releases  for  Summersville  and 
Sutton  are  100  and  75  cfs,  respectively. 

4.2.5  Rates  of  Change 

Mo  rate-of -change  criteria  were  specified  in  project  manuals  for  Sutton  or 
Summersville.  For  the  model,  these  locations  were  assigned  a  rate-of-change 
equal  to  the  maximum  channel  capacity  per  hour  at  that  location. 

4.2.6  Emergency  Pre-Release 

Whenever  forecasts  during  an  event  indicate  that  the  storage  remaining  in 
the  Kanawha  reservoirs  will  not  be  adequate  to  hold  inflows  while  continuing 
to  operate  and  keep  downstream  flows  at  channel  capacity,  larger  releases  are 
required  so  that  maximum  flood  control  storage  is  not  exceeded.  The  larger 
releases  are  minimized  by  releasing  the  excess  inflow  values  over  the  longest 
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possible  period  of  time.  The  pre-release  procedure  is  an  option  in  HEC-5Q  and 
was  included  in  the  Kanawha  model. 


4.2.7  Streamflows  and  Reservoir  Discharges 

Daily  streamflow  data,  for  inflow  to  and  discharge  from  both  reservoirs, 
were  obtained  from  the  District  office  records.  Daily  streamflows  at  most 
other  control  points  were  obtained  from  the  USGS  WATSTORK  data  system.  All 
control  points  without  observed  data  (e.g.,  London  and  Harmet  Locks  and  Dams 
for  all  periods,  and  Clay  and  Frametown  for  some  periods)  were  calculated  by 
HRC-5Q  to  minimize  negative  local  flows. 

Using  the  observed  or  calculated  flows  at  all  control  points,  the  local 
subbasin  flows  were  calculated  by  HEC-5Q.  The  model  routes  the  flow  from  an 
upstream  control  point  to  the  next  downstream  control  point  and  combines 
similarly  routed  flows  from  major  tributaries  (e.g..  Elk  or  Gauley  Rivers). 

This  routed  flow  was  subtracted  from  the  downstream  control  point  flow  to 
determine  the  local  subbasin  flow.  Some  negative  local  flows  will  occur  but 
they  are  minimized  by  adjusting  routing  criteria. 

4.2.8  Reservoir  Storages 

The  reservoir  inflows  are  routed  through  each  reservoir  with  known 
discharge  and  estimated  evaporation  to  determine  a  calculated  reservoir 
storage  and  associated  elevation.  The  calculated  elevation  was  compared  to 
the  observed  elevation.  Accuracy  of  the  results  is  discussed  in  Section  4.2.9. 

4.2.9  Evaporation  and  Precipitation 

The  net  evaporation  data  (evaporation  minus  precipitation)  used  for  HEC-5Q 
was  obtained  from  EM  1110-2-1701  [U.S.  Army,  1952]  for  a  critical  period  type 
analysis.  Although  the  data  is  not  necessarily  appropriate  for  the  actual 
historical  period,  the  simulations  demonstrate  average  daily  reservoir 
elevation  accuracy  to  within  ±0.7  feet,  on  the  average,  at  both  reservoirs 
for  1976,  1979,  and  1980. 

Demonstration  of  the  validity  of  this  type  of  data  allows  its  use  for 
real-time  prediction. 


4.2.10  Weather  Conditions 

The  meteorological  data  was  obtained  from  the  NOAA  National  Climatic 
Center  in  Asheville,  North  Carolina.  Four  stations  were  obtained  for  the 
period  1961-1983  with  some  periods  missing  at  some  stations.  The  stations 
obtained  were  Huntington  (elev.  827),  Berkley  (elev.  2504),  Charleston  (elev. 
939)  and  Elkins  (elev.  1992).  The  HBC-5Q  model  is  able  to  use  any  one  of  the 
stations  in  any  part  of  the  basin  analysis. 
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4.2.11  Intake  Level  and  Discharge  Capacity 


The  physical  description  of  the  intake  design  and  the  discharge  capacity 
were  obtained  from  the  reservoir  regulation  manuals.  Sutton  and  Sutnmersville 
Lakes  intake  data  are  sho%m  in  Table  4.1. 


Table  4.1 


Intake  Elevations  and  Discharge  Capacities 


Outlet  Type 


Sutton 

(feet/cfs) 


Low  flow  outlet 
Flood  flow  outlet 
Emergency  spillway 


910.0/  1,500 
830.0/  23,900 
972.0/222,000 


4.2.12  Reservoir  Width  and  Length 


Summersville 

(feet/cfs) 

1400.0/  318 

1398.4/  17,650 
1710.0/412,000 


The  width  and  effective  length  of  the  reservoir  were  estimated  based  on 
plan  view  diagrams  in  the  reservoir  regulation  manuals.  The  effective  lengths 
were  eventually  decreased  by  one-half  for  improved  reproduction  of  observed 
data.  Final  values  used  for  widths  and  lengths  for  Sutton  and  Summersville 
Lakes  are  shown  in  Table  4.2. 


Table  4.2 


Reservoir  Widths  and  Effective  Lengths 


Reservoir 

Eff.  Length 

Widths/Elevations 

(feet) 

(feet/feet) 

Sutton 

34,000 

275/  825 

300/  850 

325/  870 

350/  895 

375/  910 

400/  922 

425/  940 

450/  960 

475/  972 

500/  976 

525/  980 

550/  985 

575/  990 

600/1000 

625/1001.5 

650/1003 

675/1005 

700/1009 

Summersville 

21,000 

400/1391 

600/1535 

800/1555 

1000/1575 

1200/1600 

1400/1620 

1600/1652 

1800/1680 

2000/1710 

2200/1712 

2400/1714 

2600/1716 

2800/1722 

3000/1728 

3200/1735 

4.2.13  Reservoir  Clarity  and  Diffusion 

The  clarity  of  Sutton  Reservoir  is  less  than  that  of  Summersville,  and 
their  respective  visibilities  are  assumed  to  be  8  and  15  feet  for  average 
study  period  secchi  disk  readings.  These  values  are  based  on  annual  ranges 
obtained  from  the  District  of  5-15  feet  and  10-30  feet,  respectively. 

Other  required  model  data  that  effect  solar  radiation  calculations  are 
XQPCT  and  XQDEP  as  defined  in  the  HEC-5Q  Users  Manual  [HEC,  1984a].  These  are 
shown  in  Table  4.3. 
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Reservoir  diffusion  using  the  stability  method  was  initially  estimated 
from  the  HEC-5Q  Users  Manual  (page  43}  and  then  calibrated  to  reproduce 
observed  reservoir  temperature  profiles  obtained  from  the  EPA  STORET  data 
system.  The  final  values  for  diffusion  with  the  stability  method  are  also 
shovm  in  Table  4.3. 


Table  4.3 


Solar  Radiation  and  Diffusion  (Stability  Method)  Data 


Reservoir 

XQPCT 

(fraction) 

XQDEP 

(feet) 

Al 

GSWH 

A3 

Sutton 

.6 

3 

1.0x10-* 

0,9x10-6 

-.7 

Summersville 

.6 

6 

2.3x10-* 

0 . 9x10-6 

-.7 

Slight  improvements  in  reproduction  were  detected  by  changing  to  the  wind 
mixing  diffusion  method  and  using  GMIN.  GSWH,  Al,  A2  and  A3  equal  to  0, 

10~^,  10“^,  4.6  and  2  x  10"^  respectively. 

4.2.14  Decay  Coefficients 

The  decay  coefficients  used  in  the  reservoirs  and  river  network  were  the 
recommended  model  default  values.  No  attempt  was  made  to  modify  these  values 
during  the  model  calibration. 

4.2.15  Stream  Network 

The  stream  network  was  subdivided  into  stream  reaches.  A  stream  reach  is 
defined  between  each  pair  of  control  points  shown  in  Figure  2.1.  The  reach 
numbers,  the  upstream  and  downstream  control  points  numbers  and  river  mile 
identification  for  each  reach,  and  the  element  length  (i.e.,  distance  between 
computation  nodes)  are  shown  in  Table  4.4. 

Table  4.4 

Stream  Reach  Network 


Reach 

Upstream 

Downstream 

Element 

Number 

CP 

RM 

CP 

RM 

Length 

(#) 

iil 

(mile) 

Lil 

(mile) 

(miles) 

1 

954 

101.1 

966 

80.8 

1.015000 

2 

966 

80.8 

968 

52.5 

1.010714 

3 

968 

52.5 

970 

25.8 

0.988889 

4 

970 

25.8 

980 

0 

0.992308 

5 

959 

34.3 

920 

6.9 

1.014815 

6 

920 

6.9 

930 

0 

0.985714 

7 

845 

159.0 

930 

94.5 

1.954545 

8 

930 

94.5 

20 

82.8 

0.975000 

9 

20 

82.8 

10 

67.7 

1.006667 

10 

10 

67.7 

980 

54.3 

1.030769 

11 

980 

54.3 

1 

31.1 

1.008696 
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4.2.16  Cross  Sections  and  Energy  Slope 

The  cross  section  data  was  obtained  from  the  District  office.  Some 
sections  were  in  graphical  form  and  some  were  in  HEC-2  [HEC,  1982]  format  on 
magnetic  media.  The  energy  slope  for  all  river  reaches  was  estimated  from 
hydraulic  slopes  provided  on  the  navigation  river  atlas  [U.S.  Army,  1985]  or 
on  copies  of  profile  sheets  from  past  studies  in  the  District  files. 

4.2.17  Water  Quality  Concentrations 

Water  quality  concentrations  for  all  reservoir  inflows  and  initial 
reservoir  profile  conditions  were  obtained  from  the  U.S.  EPA  STORET  data 
system,  as  were  other  upstream  boundary  flows  (i.e.,  at  Hinton),  all  local 
subbasin  inflows  and  numerous  points  in  the  New,  Kanawha,  Gauley  and  Elk 
Rivers.  All  of  the  data,  except  for  the  later  points,  were  considered  to  be 
only  initial  values  which  were  ultimately  calibrated  to  final  values  by 
reproducing  the  results  along  the  stream  network.  The  adopted  concentration 
and  seasonal  variations  (but  constant  pattern  from  year-to-year)  are  shown  in 
Table  4.5.  The  HEC-5Q  model  uses  linear  interpolation  of  the  input  for 
seasonal  variations.  Demonstration  of  the  validity  of  using  the  same  seasonal 
values  for  each  year  simulated  allows  the  model  to  be  used  for  real-time 
prediction. 


Table  4.5 


Water  (Quality  Inflow  Concentrations 


Location 

Date 

Temp* 

EC 

BOD 

DO 

i!Ci 

(micromho) 

(mK/£) 

(X  Sat.) 

Sutton  Lake  Inflow 


6/1 

-4 

50 

2 

70 

9/1 

0 

100 

2 

70 

10/31 

0 

110 

2 

70 

Summersville  Lake  Inflow 

6/1 

-8 

45 

3 

90 

9/1 

-5 

80 

3 

80 

10/31 

-5 

80 

3 

80 

Hinton  Local 

6/1 

-2 

160 

2 

85 

9/1 

+3 

160 

2 

80 

10/31 

+3 

160 

2 

80 

Frametown  Local 

6/1 

-2 

130 

3 

55 

9/1 

-2 

80 

2 

55 

10/31 

-2 

80 

2 

55 

*  Temperature  departure  in 

“C  from  the 

equilibrium  temperature. 
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Table  4.5  (cent.) 


Location 

Date 

Temp* 

EC 

BOD 

DO 

(°C) 

(micromho) 

(%  Sat.) 

Clay  Local 


6/1 

-2 

20 

2 

35 

9/1 

-2 

20 

2 

35 

10/31 

-2 

20 

2 

35 

Queen  Shoals  Local 

6/1 

-2 

420 

5 

30 

9/1 

-2 

420 

5 

30 

10/31 

-2 

420 

5 

30 

Belva  Local 

6/1 

-4 

450 

10 

100 

9/1 

-4 

450 

10 

100 

10/31 

-4 

450 

10 

100 

Kanawha  Falls  Local 

6/1 

-4 

120 

4 

85 

9/1 

-1 

120 

4 

75 

10/31 

0 

120 

4 

75 

London  L/D  Local 

6/1 

-4 

400 

4 

75 

9/1 

-1 

400 

4 

75 

10/31 

0 

400 

4 

75 

Marmet  L/D  Local 

6/1 

-4 

400 

4 

95 

9/1 

-1 

400 

4 

95 

10/31 

0 

400 

4 

95 

Charleston  Local 

6/1 

-2 

800 

50 

20 

9/1 

-1 

800 

50 

20 

10/31 

0 

800 

50 

20 

Winfield  L/D  Local 

6/1 

-2 

160 

10 

90 

9/1 

+1 

160 

10 

85 

10/31 

+2 

160 

10 

85 

*  Temperature  departure  in  °C  from  the  equilibrium  temperature. 

4.2.18  Control  Point  Water  Quality  Targets 

The  water  quality  targets  at  all  control  points  were  provided  by  the 
District  based  on  locally  desired  water  quality  conditions.  Maximum 
acceptable  temperature  of  24 °C  and  minimum  acceptable  dissolved  oxygen  of 
6  mg/&  were  used  as  target  values  at  all  control  points  except  Sutton, 
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Frametown,  Clay  and  Queen  Shoals.  These  control  points  on  the  Elk  River  were 
control  for  minimizing  acceptable  temperatures  of  24®C  which  causes  the  use  of 
the  high  rise  outlet  at  Sutton.  These  values  are  used  when  making  operational 
executions  which  allow  the  program  to  decide  gates  settings  and/or  reservoir 
discharges . 

4.3  DATA  MAHIPULATIOH 

Two  types  of  data  (i.e.,  channel  geometry  and  weather)  require  independent 
manipulation  with  utility  programs  prior  to  input  to  the  HEC-5Q  model.  The 
channel  cross-section  data  is  normally  found  in  either  graphical  form  or 
already  processed  into  a  format  used  for  the  HEC-2  program  [HEC,  1982].  When 
only  the  graphical  form  was  available,  a  utility  program  was  first  used  in 
conjunction  with  a  graphics  tablet  to  process  the  data  into  the  HEC-2  format. 
Once  the  cross-sections  are  available  in  the  HEC-2  format,  they  are  input  to 
program  GEDA  [HEC,  1981]  which  processes  them  for  direct  input  to  HEC-5Q. 

The  NOAA  National  Climatic  Center  weather  data  was  obtained  in  their  C0144 
format.  The  file  was  processed  through  a  utility  program  WEATHER  [HEC,  1986]. 
The  output  file  from  WEATHER  is  input  to  the  HEATX  utility  program  [U.S.  Army, 
1977]  which  processes  the  data  for  direct  input  to  HEC-5Q. 

Procedures  for  executing  GEDA,  WEATHER  and  HEATX  are  documented  in  Appendix 

0. 

4.4  MODEL  CALIBRATION  CONSIDERATIONS 

The  objective  of  calibration  was  to  make  the  model  simulate  the  actual 
regulation  of  the  Kanawha  Basin  system  as  closely  as  possible. 

The  calibration  process  generally  consisted  of  assigning  estimated  values 
to  unmeasurable  or  unmeasured  data.  If  the  results  of  the  simulation  were  not 
appropriate  or  satisfactory,  other  values  were  assigned  and  the  simulation  and 
evaluation  were  repeated.  In  circumstances  where  the  relationships  between 
variables  and  results  were  ambiguous  due  to  interactions,  interacting  variables 
were  calibrated  on  separate  steps.  Calibration  on  a  full  range  of  possible 
events  was  not  possible  because  of  the  limited  availability  of  calibration 
data.  Instead,  calibration  focused  upon  three  of  the  four  events  chosen  by  the 
Huntington  District.  The  fourth  event  was  used  for  verification  data.  All 
four  events  were  representative  of  the  more  frequent  low  flow  condition. 

4.4.1  Data 

Data  for  October  1976  and  June  through  September  1979  and  1980  were  used 
for  calibration  purposes.  Data  from  the  events  reflected  typical  low  flow 
reservoir  regulation  periods  very  similar  to  the  type  of  periods  during  which 
the  model  will  be  used  as  a  forecasting  tool. 

Data  for  the  calibration  events  were  average  daily  flows  and  instantaneous 
water  quality  samples.  Weather  data  was  daily  average  except  cloud  cover  which 
was  averaged  over  the  daylight  hours  only  (this  parameter  is  used  for 
attenuation  of  solar  radiation) . 


4.4.2  Water  Quantity  Calibration 

The  water  quantity  aspects  of  calibration  effect  both  the  reservoir 
continuity  and  the  routing  and  combining  of  discharges  front  local  subbasin 
flows.  The  water  quantity  part  of  the  data,  with  the  exception  of  evaporation, 
had  been  previously  calibrated  and  documented  in  the  Real-Time  Flood 
Forecasting  Report  [HRC,  1984b]. 

4. 4. 2.1  Evaporation 

The  evaporation  data,  used  in  the  real-time  model,  was  modified  (because  a 
low  flow  season  has  more  evaporation  than  flood  events)  by  values  calculated 
according  to  EM  1110-2-1701  lU.S.  Army,  1952],  and  the  reservoir  continuity  was 
checked  by  comparing  reservoir  pool  elevations  during  all  three  study  periods. 
The  results  were  accurate  to  within  ±0.65  feet  (average  for  all  three 
events) . 


4 . 4 . 2 . 2  Channel  Routing 

No  routing  is  necessary  for  daily  flows  over  reaches  having  travel  times 
which  are  much  less  than  one-day.  Only  the  one  reach  from  Hinton  to  Kanawha 
Falls  had  a  travel  time  exceeding  a  half  day  (i.e.,  about  18  hours)  and  routing 
criteria  were  estimated  using  reproduction  of  observed  streamflow  data  as  the 
evaluation  criteria. 

4.4.3  Water  (Quality  Calibration 

The  water  quality  calibration  involves  adjustment  of  four  reservoir 
diffusion  coefficients  (Al,  A2,  A3  and  GSWH)  and  three  solar  energy 
distribution  coefficients  (EDMAK,  XQPCT,  and  XQDEP) .  Sensitivity  studies  were 
used  to  select  the  most  appropriate  weather  station  to  represent  the  weather  at 
each  reservoir  and  within  each  stream  reach.  Since  local  subbasin  inflow 
quality  is  relatively  unknown,  sensitivity  studies  were  also  used  to  adjust 
those  values.  All  of  the  above  are  adjusted  through  numerous  trials  to  best 
reproduce  observed  reservoir  profiles  and  observed  time  series  data  at  river 
sampling  locations. 

4. 4. 3.1  Diffusion  Effects 

The  GSWH  and  A3  are  minimum  and  maximum  diffusion  coefficients.  The  higher 
the  Al  value,  the  more  heat  that  is  transferred  up  to  a  maximum  value  of  A3. 

The  A2  value  impacts  the  shape  of  the  diffusion  function  with  depth.  The  A2 
value  is  recommended  to  be  held  constant  at  4.6.  Magnitudes  of  GSWH,  Al,  and 
A3  derived  by  successive  trials  for  both  Sutton  and  Summersville  are  10~^, 

10"^,  2  X  10“^  respectively. 

4. 4. 3. 2  Solar  Radiation  Distribution 

The  EDMAK  (Secchi  disk)  controls  the  depth  of  light  penetration  which 
effects  the  distribution  of  the  solar  radiation  which  then  penetrates  the  water 
surface.  The  larger  the  EDMAK,  the  deeper  the  solar  energy  penetration. 

Secchi  disk  readings  are  field  observations  of  reservoir  clarity  performed  by 
gradually  lowering  a  black  and  vrhite  disk  to  a  depth  where  it  can  no  longer  be 
seen.  The  depth  where  it  disappears  from  sight  is  recorded. 
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The  fraction  of  the  total  solar  energy  penetrating  the  water  surface, 
which  is  absorbed  near  the  surface  of  the  reservoir,  is  defined  as  XQPCT. 

The  larger  the  value,  the  wanner  the  water  surface. 

The  depth  of  water  that  absorbs  XQPCT  is  defined  as  XQOEP.  The  deeper 
XQDEP,  the  smaller  the  actual  temperature  increase  of  the  surface  water. 

The  value  of  EDMAX  depends  on  observed  data  in  each  reservoir.  At 
Sutton,  the  EONAX  was  estimated  to  be  8  feet  and  at  Summersville,  it  was  15 
feet.  The  magnitudes  of  XQPCT  and  XQDEP  are  not  usually  measured,  but 
estimated  values  of  0.6  and  75%  of  the  top  layer  depth  were  used. 

4. 4. 3. 3  Weather  Station  Selection 

To  decide  on  the  most  representive  weather  stations  to  use,  a  sensitivity 
study  was  performed.  Initial  trials  at  each  reservoir  used  the  closest 
weather  station  to  their  geographical  location.  Similarly,  each  river  reach 
was  assigned  the  closest  weather  station.  Eventually,  the  Elkins  weather 
station  was  selected  because  it  provided  the  best  representation  of  weather 
conditions  for  both  reservoirs  and  all  the  river  reaches.  Use  of  the  Beckley 
data  provides  only  minor  variations. 

4. 4. 3. 4  Subbasin  Inflow  Quality 

The  data  from  the  STORET  system  is  the  best  starting  place  for 
determination  of  subbasin  inflow  quality.  Unfortunately,  the  data  is  sparce 
in  time  and  space.  The  data  is  normally  measured  about  once  per  month  on  a 
few  tributaries.  Even  this  limited  data  may  not  be  available  near  the 
downstream  end  of  the  tributary.  Due  to  these  limitations,  STORET  data  only 
provides  a  ‘'ballpark"  concentration  level  and  a  trend  for  its  seasonal 
variation.  Final  input  data  must  be  determined  through  calibration  or 
sensitivity  studies.  Section  4.2  of  this  report  previously  defined  the 
subbasin  inflow  concentrations  finally  adopted.  It  is  important  to  note  that 
since  they  were  kept  constant  for  all  three  years  at  all  subbasins,  the  same 
values  should  also  be  acceptable  for  real-time  prediction  use. 

4.5  CALIBRATION  RESULTS 

Using  the  calibration  coefficients  and  parameters  previously  discussed, 
the  model  is  capable  of  satisfactory  reproduction  of  the  1976,  1979,  and  1980 
observed  reservoir  and  stream  quality  data.  Evaluation  of  the  final  results 
must  be  examined  in  terms  of  all  study  years,  all  seasonal  variations  and  all 
geographical  locations.  Changing  a  coefficient  or  parameter  at  one  point  in 
time  and  space  can  seriously  effect  other  locations  and  points  in  time  (i.e., 
seasonal  variations). 

Several  examples  of  the  graphical  reproductions  are  included  in  this 
report  as  Appendix  A.  The  Sutton  (pages  A1-A8)  and  Summersville  (pages 
A9-A16)  reservoir  profiles  for  temperature  and  dissolved  oxygen  (DO)  have 
been  included  for  the  four  low-flow  months  of  either  1979  or  1980,  depending 
on  the  amount  of  observed  data  available.  The  Sutton  temperature  profile 
reproductions  are  generally  of  comparable  shape  and  within  ±2‘’C.  The 
Sutton  DO  profile  reproductions  are  of  similar  shape  except  at  the  deepest 
levels.  In  general  the  accuracy  is  with  ±1  mg/ft. . 
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Summersville  temperature  and  DO  profile  reproductions  are  all  of  similar 
shape  except  for  the  deepest  levels  and  within  similar  accuracy  to  Sutton. 

At  the  deeper  levels,  the  computed  temperatures  are  warmer  and  the  computed 
DO  higher.  The  computed  DO  cannot  reproduce  the  observed  DO  spike  which 
occurs  at  approximately  220  feet  of  depth. 

The  accuracy  of  profile  reproduction  at  both  reservoirs  is  sufficient  to 
provide  accurate  predicted  discharge  water  quality  for  downstream  analysis. 

The  Gauley  River  temperature  (pages  A17-A18)  and  DO  (pages  A19-A20) 
low-flow  time  series  for  1979  have  been  shown  for  below  Summersville  Dam  and 
near  Belva.  The  accuracy  of  reproduction  is  within  ±2°C  for  temperature 
and  ±.5  mg/9,  for  DO.  The  graphs  on  pages  A17  and  A19  support  the 
previous  statement  regarding  the  adequacy  of  the  calculated  reservoir  profile 
accuracy  in  Summersville. 

The  Elk  River  temperature  (pages  A21-A29)  and  DO  (pages  A30-A38)  time 
series  for  1980  are  shown  for  four  locations  between  Sutton  Dam  and  Frametown 
( EiM98-RM90 } ,  for  one  location  between  Frametown  and  Clay  (RM62),  for  one 
,  between  Clay  and  <}ueen  Shoals  (RM52)  and  for  three  locations  between  Queen 

i  Shoals  and  Charleston  (RM2S-RM4).  The  accuracy  of  reproduction  is  generally 

within  ±2  C  for  teRq>erature  and  ±1  mg/ ft.  for  DO.  The  graphs  on  pages 
A21  and  A30  support  the  previous  statement  regarding  the  adequacy  of  the 
calculated  reservoir  profile  accuracy  in  Sutton. 

I 

The  Hew/Kanawha  River  temperature  (pages  A39-A42)  and  DO  (pages  A43-A46} 
time  series  for  1979  are  shorn  for  four  locations  including  Hinton  (RM161), 
above  Kanawha  Falls  (RM98),  and  the  Marmet  pool  (RM73)  and  in  the  Winfield 
pool  below  Charleston  (RM46}.  The  accuracy  of  reproduction  is  generally 
within  ±2'*C  for  temperature  and  ±1  mg/ft  for  DO.  The  graphs  on  pages 
A39  and  A43  show  that  the  boundary  conditon  at  Hinton,  which  is  analyzed  as  a 
"dummy**  reservoir,  is  a  sufficiently  adequate  representation  of  the  physical 
condition. 

Graphs  for  other  years  having  iess  observed  data,  and  for  electrical 
conductance  (EC)  and  biochemical  oxygen  demand  (BOD)  have  been  furnished  to 
Huntington  District  under  separate  cover.  Because  EC  and  BOD  are  of  less 
concern  to  the  District  and  the  impracticality  of  including  about  five  times 
as  many  graphs  as  have  been  shown,  these  other  graphs  have  been  excluded.  It 
should  be  noted,  however,  that  data  for  other  time  periods  are  of  similar 
accuracy  to  those  shown  and  that  the  accuracy  of  EC  and  BOD  reproductions  are 
generally  within  ±50  micromhos  and  ±1  mg/ft  respectively. 

In  general,  the  reservoir  reproductions  are  good  and  the  river 
reproductions  are  very  good  overall.  The  calibrations  shown  in  Appendix  A 
were  accomplished  with  more  attention  paid  to  determining  inflow 
concentrations  that  can  be  used  for  future  real-time  forecasting  than 
obtaining  an  exact  reproduction  of  a  given  value  on  a  specific  date  or 
location. 

The  error  that  can  be  observed  in  the  results  shown  in  Appendix  A  is  the 
same  magnitude  of  error  expected  in  real-time  work.  As  more  experience  with 
the  model  is  gained,  additional  fine  tuning  may  result  in  small  improvements 
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in  calibration.  Any  time  the  user  decides  to  try  a  model  modification,  these 
calibration  data  should  be  rerun  and  the  new  graphical  output  compared  to 
these  results. 

4.6  MODEL  VERIFICATION 

Model  verification  is  the  best  available  proof-of-model-validity  for  use 
in  real-time  prediction  or  for  use  under  proposed  future  system  structural 
changes.  The  verification  test  must  use  all  the  calibration  coefficients  and 
parameters  as  constants  (i.e.,  these  values  cannot  be  changed  on  the 
verification  case). 

The  verification  case  chosen  by  the  District  was  June  through  September 
1983.  The  input  data  for  this  case  included  the  same  physical  description 
data  for  the  reservoirs  and  stream  channel  as  used  previously  for  the 
calibration  cases.  The  weather  data  and  observed  water  quality  data  from  the 
STORET  system  were  obtained  for  the  1983  period.  The  local  subbasin  flow 
quantity  was  determined  by  the  model  and  the  flow  quality  was  held  constant 
at  the  calibration  concentrations. 

The  inflow  to  both  reservoirs  and  the  discharges  from  them  were  obtained 
from  Huntington  District. 

A  sample  of  the  graphical  results  of  the  verification  test  case  are  shoim 
in  Appendix  B.  The  results  were  quite  satisfactory  in  the  river  network.  No 
observed  reservoir  profiles  were  available,  but  the  comparisions  of  magnitude 
and  shape  to  previous  years'  profiles  were  satisfactory.  The  order  of  the 
graphs  and  the  locations  provided  are  generally  similar  to  the  calibration 
graphs  in  Appendix  A  except  that  the  reservoir  profiles  are  not  included  and 
observed  values  were  available  at  fewer  locations. 

The  stream  temperatures  were  slightly  low  and  the  DO  slightly  high.  This 
is  a  normal  relationship  to  be  expected  due  to  the  dissolved  oxygen 
saturation  functions.  Again,  as  discussed  above,  further  fine-tuning  of  the 
model  may  be  necessary  as  more  conditions  are  examined;  but  for  prediction  of 
these  1983  verification  results,  the  reproduction  is  reasonable. 

The  model  having  passed  the  verification  test  is  available  and 
recommended  for  field  use.  Caution  must  still  be  applied  to  interpreting  the 
model  results,  as  with  any  new  product.  If  any  anomalies  are  observed, 
additional  calibration  may  be  necessary  to  fine-tune  the  model. 

4.7  OPERATIONAL  CAPABILITY 

The  data  files  for  the  four  study  years  have  also  been  used  in  an 
operational  mode  where  the  model  decides  the  reservoirs'  discharges.  As 
shown  in  Appendix  C,  the  low-flow  (regulated)  discharges  determined  by  the 
model  (dashed  line)  are  almost  identical  reproductions  of  the  observed 
discharges  (solid  line)  at  the  dam  site  and  at  downstream  control  points  for 
all  three  study  years  (1976,  1979  and  1980).  Comparison  of  calculated  water 
quality  to  observed  water  quality  is  similar  in  accuracy  to  the  calibration 
results  for  the  same  periods  and  generally  the  same  locations  provided  in 
Appendix  A. 
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The  modifications  to  the  calibration  data  files  to  develop  operational 
capability  include  removal  of  the  specified  reservoir  releases  and 
substituting  local  flows  for  observed  USGS  data  on  the  IN  records.  The  J3.6 
index  is  changed  to  -1  to  notify  the  model  of  this  change  in  input.  Other 
modifications  include  changing  the  J9.4  index  to  0  and  removing  the  G1  and  G2 
data. 

4.8  REAL-TIME  CAPABILITY 

For  forecasting  at  a  given  time  using  the  operational  capability  type 
data  file,  update  the  RL.2  storage  value  to  the  current  data,  as  well  as  all 
of  the  local  flows  on  the  IN  records  (or  ZR  records  using  OSS).  The  number 
of  days  to  be  forecasted  and  the  current  date  must  be  updated  on  BF.2,  BF.5, 
JA.2  and  JA.3.  The  reservoir  water  quality  initial  condition  profiles  (L9 
through  C7)  and  the  weather  data  (ET)  must  also  be  updated. 

If  the  IN  and  ET  records  are  defined  for  dry,  wet,  and  average  conditions 
for  the  whole  dry  season,  the  updating  of  these  records  for  each  weekly 
forecast  can  be  eliminated. 

For  forecasting  water  quality  conditions  with  specified  reservoir 
releases,  add  QA  records  after  the  IN  records  which  define  the  reservoir 
releases  for  the  forecast  period. 


5.  SUMMARY 

The  purpose  of  the  Kanawha  River  Basin  water  quality  study  was  to  provide 
a  computer  model  capable  of  evaluating  water  quality  conditions  that  occur 
throughout  the  basin  during  June  through  September  for  a  given  operation  of 
Sutton  and  Summersville  Dams. 

Within  acceptable  tolerances,  the  study  has  accomplished  this  purpose,  as 
well  as  providing  a  tool  for  better  understanding  and  management  of  the  data 
being  collected  by  various  agencies.  Better  input  boundary  conditions  (i.e., 
inflow  quality  to  both  reservoirs  and  water  quality  at  Hinton)  are  needed  for 
improved  reproductions  in  both  reservoirs  and  at  all  locations  in  the  upper 
stream  reaches  of  the  New,  Gauley  and  Elk  Rivers.  More  water  quality  data  on 
subbasin  inflows  would  be  beneficial  in  the  Clay  to  Charleston  area  of  the 
Elk  River  Basin.  In  general,  data  needs  exist  for  all  inflowing 
concentrations.  This  data  is  generally  unavailable  and  yet  crucial  to  the 
study. 

The  study  helped  define  better  procedures  for  using  and  interfacing 
various  HEC  models  and  utility  programs.  It  is  the  intent  of  HEC  to  provide 
this  report  not  only  as  a  final  project  document  but  as  a  procedure  guide  to 
other  water  quality  studies  having  reservoir  system  operations  of 
multipurpose  projects.  A  companion  reference  procedure  guide  for  HEC-5Q 
studies  is  available  from  HEC  [Willey,  1985b]. 
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ilJMl€RSV]LLE  17JUL79*-CieS  TEMP 


2EMft?86  18  18  SB 


All 


2»«!S6  18  18 


A12 


SUMhCRSVILLE  2ISEP7Qk-Q0MP  TEfl^ 


2»ff?86  18  18  3B 


A13 


2SMRR86  18  18 


.(3LH1CRSV1LLE  17JIJL7£M-Ce3  CONS. 


as  81  81  98»:l«; 


2»«»86  18  18  SB 


A16 


SUM1€RSV1LL£  21SEP79*a]MP  0X\13EN 


2»fR8ir-  19  00 


'GftJLE 


2a«<86  19  Q0 


A18 


CONP+TB'IP  CONST  ITLIEWT 


9D  61  39a*«3 


A21 


7.548  CUMP+TEMP  C0M3TITIJENT 


2S»1F«d6  19  Q6  21 


'X^U';-RM'57  0B3  TEHF 
-  96.S53  CXiMP+TEMP  C0M3TITUENT 


12  90  61  99ifc»«2 


A23 


)k(ELK-RM92  OES  TEHF 
_  92.473  COMP+TEMP  COfJSTITUEMT 


2SM=«86  19  06  21 


A2A 


>(ELKH?M30  OBS  TEMP 
_  90.443  CfJMP+lEMP  COf'GTITUENT 


A25 


90  61 


A26 


)«(ELK-RM52  083  TEMP 
_  52.006  CG8  CCHP^-TEMP  CDNST  J  TLEh/T 


90  61 


)x(ELKH^M25  0B3  TEHP 
_  25.  304  C70  CCMP^TEMP  CONST  ]TLE^^^ 


>(ELK-RM21  0B3  TEMP 
_  21.335  COMP+TEMP  CCifeilTUENT 


19  B6  21 


4.4&S  COMP+TEMP  CONST ITUEHT 


IZ  90  61  99iH€Z 


96.533  COMP+CKYGEN  COreTITUEtn- 


90  61  98t»«3 


)^LK-RM92  OES  DO 
_  92.473  COMP+OKYGEN  OOhET  [  TUErvfT 


X;elk-^-m90  oes  do 
_  90.443  COMP+n>;YGEN  r]:iM3T[TIJEMr 


2aW86  19  06 


A34 


62.102  COMP+CKYGEH  COMSTITUEinr 


KflNflWHfl 


52.036  cm  CCMP^-OXYCtN  CONST  ITI^WF 


2SM^86  19  06 


A36 


VELK-RM25  OBS  DO 

y  \ 

-  25.334  C70  CCMP^-OXYCCN  CONSTITUENT 


2»«?86  19  06 


VELK-RM21  OBS  DO 
/  \ 

-  21.335  COMP+OXYGEN  CDf-tSTlTUEflT 


29«?86  19  06 


iXl  £90 


\^EW-RM16)  CeS  TEMP 


2SM=R86  18  29  10 


A40 


'.432  COMP+TEMP  COhETITUENT 


2»f«8e.  18  29 


A42 


72b  COMP+TEMP  COf^6TITUENT 


2a-«?86  18  29  10 


),<9-JB'l-RMl6]  CBS  DO 

-  158.023  C45  COMP+OXVGEM  CONSTITUENT 


S^lRRSe  18  29  10 


97.432  COMP+CKYGEN  CDf’eTlTUEf^lT 


2EMFH86  18  29 


iMP^^YGEN  CCifE;T[TUEf-n' 


2»ffl86  18  29  10 


A46 


^FW-RM46  OBS  DO 
-  45.726  COMP+OXYGEN  0:iM3T(  TUErn- 


2BHFfi86  13  07 


CCMPfTEMP  CONST 


2BMnR86  13  07  2 


COHP+TEMP  CONSTITUENT 


2EMRR86  13  07  2 


VOXYGEN  CONST ITLEWr 


28MRR86  13  07 


COMP-fOXYGEN  CCNSTirUENT 


2BM^86  13  07 


B6 


'  JUNaa  I  JUL83  I  RUG83  I  SERBS  1  0CTS3 

)h(ELK-RM25  OES  temp 
-  25.304  C70  CCMPfTEMP  CONST  ITLElfT 


2a-fR86  13  OT  2 


B7 


N^^ELK-RMi  DBS  TEMP 
-4.4GB  CONP+TDIP  CONSTITUENT 


2eMPR86  13  D1  2 


),(ELK-RM98  OBS  DO 

-  98.962  COMP+OXYGEN  OO^GT[TUE^^• 


2EMRR86  13  07 


B9 


)»(ELK-RM25  DBS  00 

_  25.304  C70  CGMP*-0XYGEN  CONST  ITLENF 


2fiMRR8e  13  07 


^LK-RM4  OBS  DO 
_  4.4GB  COMPfOX'i'GEN  CCTIST I TUENT 


28MFfi86  13  QT  2 


:5  COMP+lEt'P  CONSimJEhJT 


31^flR86  10  S5 


97.432  CCiMP+TEMP  COMSTITUENT 


31MqR86  10  SS  09 


B13 


SUfflSC  10  55  09 


B15 


^NFf'l-RMSe  OB3  TEMP 
_  S8.503  C)0  T10  C0NP*TD’1P  CONSTITUENl 


2a‘fiR86  13  07 


5  OOMP+OX'YGEN  CONSTITUENT 


3UPR36  10  55  09 


83.287  C20  T  9  CONP-.OWGEW  CCtlST  I  WEN! 


3UfflS6  10  55  B9 


B21 


?..233  C)0  T10  COMP-fOVi'GEN  CCtgSTI TUENT 


31Kffl86  10  55 


B22 


APPBiDIZ  C 


Q8APHICAL  RESULTS  OP  OPERATIONAL  CAPABILITY 


TABLE  OF  CONTENTS 


Location 

Parameter 

Page 

Sutton 

Reservoir  Release 

Cl 

SiaaMrsville 

Reservoir  Release 

C4 

Gaulop  River:  Belva 

Regulated  Flow 

C7 

Elk  River: 

Fraowtoun 

Regulated  Flow 

CIO 

Clay 

Regulated  Flow 

C13 

Queen  Shoals 

Regulated  Flow 

C14 

New/Ranawha  River: 

KanaWha  Falls 

Regulated  Flow 

C17 

Charleston 

Regulated  Flow 

C20 

Sutton 

Tes^erature 

C23 

Dissolved  Oxygen 

C27 

SuBsnersville 

Temperature 

C3I 

Dissolved  Oxygen 

C35 

Gauley  River:  RII344  -  RMS 

Temperature 

C39 

Dissolved  Oxygen 

C41 

Elk  River:  RUBS  -  RM4 

Tesverature 

C43 

Dissolved  Oxygen 

C52 

New/Kanauha  River:  BM158  -  BM46 

Temperature 

C61 

Dissolved  Oxygen 

C65 

Cl 


C2 


C3 


C5 


ILLC 


ri.iMh'  FMiU-Prr, 


C9 


CIO 


ffrm'town  aiff  fL.MW-F.'F.r. 


Cll 


C12 


C13 


C14 


C15 


OLED'J  'i.HijHL:.  COMP  FLOl'*-REG 


C16 


KhWrtiHH  fhlls  roMp  n.Li/-Pfj;- 


C18 


H^FlF.' i  Hfj  rrnp  noii-rtG 


C21 


C22 


2BRPR86  15  09  32 


C23 


2BnPR86  15  09  32 


C24 


Sl/TTON  08JUL80+COHP  TEMP 


2BFPR86  15  09 


C25 


SUTTON  04flUG80+COHP  Te-IP 


2BfiPR86  15  09 


2BfiPR86  15  09 


04JUN80+CONP  OXrGEN 


2BfiPR86  15  09  32 


C28 


.  Sl/TTON  08JULe0+CONP  O)0i'GEN 


2BPPR86  15  09  32 


C29 


SUTTDN  04RUG80-t-CONP  OXYGEN 


20fiPR86  15  09  32 


C30 


SLITTON  23SEP80+CONP  OXYGEN 


tl  ft  SI 


C3l 


v'Sl.HfCF.-'jV )  LLE  14  .lUN7:;n-Ce' 


15  31  11 


SI>11€F'Vm'1LLE  16f(UG7'4«.iD3  TEMP 


U  IE  SI 


C34 


SLMr-ERSVILLE  21SEP79^mMP  lEIP 


C35 


C36 


Cl  ct  SI 


C37 


SLMf€R'iV]LL£  lt:fllJi:.7'3^a:iMF'  iVrfTV.DC 


C38 


SLnftRSVlLir  21SEP794-03MP  Cixvir-D'^ 


2EPFfl8e.  12  52  SB 


C39 


I  JUNig  I  JUL79  I  flU3?3  I  SEF79  I  OCTTh 

''.(GFULEY-RHS'l  CBS  TEMP 
_  33.733  C59  CCMPfTEMP  CONST  ITLEf'JT 


2EFPR36  12  52  Sf 


;:.422  CONP+TB'IP  CONSTITUENT 


2£fPR86  J2  52 


/ijflUL£Y-RH34  CBS  DO 


2B¥«36  12  52  3E 


C42 


8.42  CONP-tOCTGEN  CCMSTIFUENT 


28FrR86  14  00  3 


C43 


)>(ELK-RM98  OBS  TEHP 
_  97.548  COMP+TEMP  CONBTITUEKT 


2BFrR86  14  00  3 


)k(ELK-RM97  OBS  TEHP 
-  9G.B33  COMP+TEMP  CONSTITUENT 


2BPPR86  14  00  3 


^LK-1?M90  0B5  TEMP 
_  90.443  COMP+TEMP  CONSTITUENT 


KflNflWHfl 


C47 


)k(ELK-RM62  C  BS  TEHP 
_  62. 102  COMP+TEMP  CONSTITUENT 


2BPPR86  14  00  3 


C50 


XELK-RM21  08S  TEMP 

_  21.335  COMP+TEMP  COieTITUENT 


2em86  14  00  34 


C51 


)h(ELK-«M4  DBS  TEMP 
-  4.465  CONP-^TEMP  CONSTITUENT 


2efiPR86  14  00  3 


)<^LK-»?M98  OBS  DO 
-  97.548  COMP+OXYGEN  CONST  [TUEf^fT 


2KrRe6  14  Q0  3 


C53 


OBS  D0 

_  9G.533  COMP+OXYGEN  C01GTlTUEf>fT 


2BFrR86  14  00  3 


)h(ELK-RM92  OBS  do 

_  92.473  COMP+CKYGEN  CO^BT(TUE^^n• 


2BRPR86  14  00  3 


C55 


)k;^LK-RM90  OBS  DO 
_  90.443  COMP+CKYGEM  ODNSTITUENT 


2BflPR86  14  00  3 


),(EU<-RM62  OBS  DO 
-  62.102  COMP+CKYGEM  CONST (TUEMT 


2BFPR86  14  00  3 


C58 


)h(ELK-1?M25  OBS  do 

_  25.304  C70  CCMP*-OXYGEN  CONST  ITLEWF 


2BRPR86  14  00  3i 


C59 


)<(ELK-^M21  OBS  DO 
_  21.335  COMP+OKYGEN  CDNSTITUEMT 


hwhh 


i.’Gi 


H.HNHWflH 


,•  I IK-.  rpHP 


C65 


if-Gi'ii  ]  i 'iJ-Ji;:!  j> .jwii.:! 


ij-rjfii  1  IS  4i;o  tG';iA '  jj  . 


C68 


APPEITDIX  D 


UTILITY  PROGRAM  PROCEDURES 


APPEIIDIX  D 


Utility  Program  Procedures, 


GEDA  EXECUTION  PROCEDURE 


1.  Obtain  HEC-2  type  cross-sectional  data  files  either  from  past  studies  or 
if  in  graphical  form  convert  to  a  magnetic  file  using  a  utility  program  and  a 
graphics  tablet  (available  from  HEC  but  is  currently  not  in  a  generalized, 
documented  form) . 

2.  Prepare  input  file  as  shown  in  GEDA  Users  Manual  [HEC,  1981].  The  bulk  of 
the  input  file  will  be  the  HEC-2  type  cross-sectional  data  but  be  careful  to 
edit  the  HEC-2  data  card  types,  according  to  the  guidance  in  the  GEDA  manual. 

3.  To  execute  GEDA,  type: 

GEDAX  INPUT=filename.IN  OUTPUT=f ilename . OUT  PUNCH=f ilename.PU 
(Note:  if  necessary  for  space  use  $M0  PS=500) 

4.  The  "f ilename.PU"  from  GEDA  needs  some  editing  before  it  is  ready  to  merge 
with  HEC-5Q  input.  The  following  steps  must  be  performed: 

(a)  change  the  record  identification  from  A3  to  S3  in  columns  1  and  2. 

(b)  The  second  field  of  numbers  should  be  deleted. 

(c)  The  first  field  of  numbers  should  be  changed  to  river  mile 
identification  if  the  HEC-2  type  input  used  identification  not  related  to 
river  mile.  This  step  can  be  accomplished  easily  by  a  utility  program.  The 
first  section  number  should  be  the  most  upstream  cross-section  river  mile. 

Each  cross-section  river  mile  thereafter  should  be  reduced  by  the  element 
length  (or  the  reach  length  divided  by  the  number  of  elements  in  the  GEDA 
run) .  The  last  cross-section  number  should  be  exactly  the  river  mile  of  the 
downstream  end  of  the  reach. 

(d)  Control  point  numbers  should  be  inserted  into  field  1,  making  the 
cross-section  numbers  move  to  field  2,  etc.  Each  cross-section  will  have 
either  its  own  control  point  number  or  the  number  of  the  next  downstream 
control  point  on  the  first  record  of  each  section. 

(e)  If  stage  control  (i.e.,  use  of  a  rating  curve)  in  HEC-5Q  is  important 
at  any  given  cross-section,  field  8  of  the  S3  record  should  contain  the  flow 
rate  corresponding  to  the  elevation  in  field  3.  If  this  field  is  blank, 
normal  depth  is  assumed  unless  a  flat  surface  projected  upstream  from  a 
downstream  section  using  stage  control  provides  a  larger  depth. 

5.  All  of  the  tasks  in  Step  4  can  be  easily  accomplished  by  use  of  a  utility 
routine  (available  from  HEC  but  is  currently  not  in  a  generalized,  documented 
form) . 
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WEATHER  EXECUTION  PROCEDURE 


1.  Physically  load  magnetic  tape  onto  mainframe  tape  drive. 

2.  Type  on  terminal: 

RS  14=TAPE  WAIT 
BLOCK 

TAPE  TO  DISK 

Account  Number*filename. IN  80  80  12 

(Wait  for  delay  while  tape  is  being  read) 

$EOF 
FR  14 

3.  "Filename. IN"  now  contains  all  the  years  of  data  that  were  on  the  original 
magnetic  tape.  Edit  "filename. IN"  to  include  only  those  years  of  interest 
(this  step  is  optional  to  save  computer  time).  Add  a  header  record  as  shown 
in  the  WEATHER  Users  Manual  [HEC,  1986]. 

4.  To  execute  WEATHER,  type; 

*SAUF77.W  WEATHR 
$VU.BR  WEATH,  PR  PX  OW  OD 
LIB  HECN*HECLIB  *LIBERY 
BEGIN 

WEATH  INPUT=filename. IN  OUTFILE=filename.OUT 

5.  Proceed  to  HEATX  procedure  with  "filename. OUT". 


HEATX  EXECUTION  PROCEDURE 


1.  Prepare  input  file  as  shown  in  "Thermal  Simulation  of  Lakes"  Users  Manual 
lU.S.  Army,  1977]. 

"Filename. OUT"  from  WEATHER  provides  card  types  8-11  in  the  proper 
format.  Edit  "filename. OUT"  to  include  card  types  1-7.  Rename  file 
"filename. IN". 

2.  To  execute  HEATX,  type; 

HEATX  INPUT=f ilename. IN  0UTPUT=f ilename.OUT  PUNCH=f ilename . PU 

3.  The  punch  file  is  input  to  HEC-5Q.  Merge  "filename. PU"  into  the  HEC-5Q 
input  file. 
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